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Abstract: We study the two-dimensional reduction of the Michaelis-Menten reaction of 
enzyme kinetics. First, we prove the existence and uniqueness of a slow manifold between 
the horizontal and vertical isoclines. Second, we determine the concavity of all solutions in 
the first quadrant. Third, we establish the asymptotic behaviour of all solutions near the 
origin, which generally is not given by a Taylor series. Finally, we determine the asymptotic 
behaviour of the slow manifold at infinity. To this end, we show that the slow manifold can 
be constructed as a centre manifold for a fixed point at infinity. 
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1 Introduction 

The Michaelis-Menten (more accurately the Michaelis-Menten-Henri) mechanism is the sim- 
plest chemical network which models the formation of a product through an enzymatic catal- 
ysis of a substrate. See, for example, P, [HI [13, [25], Chapter 1 of [14J, and Chapter 10 of 
[16] . In particular, an enzyme reacts with the substrate and reversibly forms an intermediate 
complex, which then decays into the product and original enzyme. Symbolically, 

S + E ^ C ^ P + E, 

k-i 

where S stands for substrate, E stands for enzyme, C stands for complex, and P stands 
for product. By the Law of Mass Action we have a set of four differential equations for the 
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concentrations s, e, c, and p\ 



s = k-ic — kise, 

e = + k2)c — kise, 

c = kise — (A;_i + k2)c, 

p = k2C. 

Taking initial conditions s(0) = sq, e(0) = eo, c(0) = 0, and p(0) = 0, we have e = eo — c and 
two independent differential equations: 

s = k^ic — kis{eo — c), 

c = A;is(eo - c) - {k_i + k2)c. 

To simplify even further, the Quasi-Steady-State Approximation (QSSA) was introduced 
by Briggs and Haldane[3]. This takes c = to hold after a short time, giving 

CqS , . . k2eos 

and —s = p 



where '■= is called the Michaelis constant. This is generally thought to be valid 

when cq <C Sq. The expression for p gives a measure of the velocity of the reaction. Experi- 
ments have been fitted to the quasi-steady-state approximation. This type of approximation 
is often used to simplify other chemical kinetics systems including those of different and more 
complicated enzyme reactions which may involve inhibition or cooperativity effects [5| [H| [T5] . 

Another approximation, though less common than QSSA, is the (rapid) Equilibrium 
Approximation (EA), which originated with Henri and was popularized by Michaelis and 
Menten. This takes s = to hold after a short time, giving 

cqs . k2eos 

and p 



Ks + s Ks + s' 

where Kg := 

The validity of QSSA can be examined from several points of view. First the equations 
are written in a non-dimensional form, introducing a small positive parameter e. This can be 
done is several different ways|71 In the traditional scaling^lQ HSl s = —. Following 
[2211231 [21], we have 

X = —X + {1 — r])y + xy, y = e~^{x — y — xy), (1) 

where x := , ^^f, is a scaled substrate concentration, -y := is a scaled complex concen- 
tration, t := A;ieoT is a scaled time, r is the original time, ' = ^, and eo := e(0) + c(0). The 
parameters are given by e := ^ ^^^^"^^ and rj := ^ ^^^^ . Note that 5 > and < ?] < 1. 

When e is small, the system ([1]) is a singular perturbation problem. A matched asymp- 
totic analysis yields the QSSA as the zeroth-order term in the outer expansion [16 j. The 
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correctness of this analysis is proved using Tikhonov-Levinson theory [20]. Exphcit bounds 
on the approximations have been obtained for small e^E\. Centre manifold theory[l] and 
geometric singular perturbation theory [T3] have been applied to give an invariant manifold 
A^£, called a slow manifold, within distance 0{e) of A4o, the quasi-steady-state manifold 
y = j^. Trajectories approach the slow manifold exponentially fast and then evolve along 
it at a slower rate. 

Several chemists have observed and theoretically investigated slow manifolds which at- 
tract other solutions. In general slow manifolds are not uniquely defined. In two-dimensional 
cases, Eraser and Roussel[8l [221 ISSl [21] take as a slow manifold the solution between the 
horizontal and vertical isoclines, which are the quasi-steady-state and the equilibrium ap- 
proximations, respectively. Roussel, for example, provided a heuristic argument based on 
antifunnel theory that there is indeed such a solution[23]. Davis and Skodje[6] take as a slow 
manifold the trajectory joining a saddle at infinity and a stable node, approaching in a slow 
direction. This paper was motivated by the work of Eraser and Roussel. 

Occasionally, we may refer to the system ([1]) in the compact form 

x = g(x). (2) 
We will also work with the one- dimensional version of ([1]), given by 

y' = fix,y), (3) 

where ' = Explicitly, 

^ \^ e'^{x — y — xy) 

In this paper, we do not need to assume that e is small. The focus is on the behaviour 
of solutions in the phase plane, that is, considering y as a function of x. In ^ we give 
the basic phase portrait of ([T]) in the first quadrant and the linearization at the origin. In 
§31 we describe the isocline structure which is exploited in subsequent sections. In §11 we 
prove the existence and uniqueness of the slow manifold, which we denote by A^, between the 
horizontal and vertical isoclines. These were discussed in a more informal way by Eraser (see, 
for example, [8]). In §5l we determine the concavity of all solutions except the slow manifold 
by analyzing an auxiliary function. In §^1 we determine the behaviour of solutions near the 
origin by using Poincare's Theorem (see, for example, [23 p. 190); one-dimensional solutions 
y{x) are generally not given by a Taylor series, which has sometimes been assumed. This 
analysis applies to any two-dimensional system with a Hurwitz-stable equilibrium point. In 
§71 we determine when solutions enter Ei, which is a region bounded below by the horizontal 
isocline and above by the isocline for the slope of the slow manifold at the origin. In §Hl we 
establish properties of the slow manifold: concavity, monotonicity, and asymptotic behaviour 
at the origin and infinity. Einally, in §9] we state some open questions. 



and f{x,y):-- 



X — y — xy 



e[—x + (1 — rj)y + xy] 
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Figure 1: A phase portrait for ([T]) for £ = 5.0 and 77 = 0.8. 



2 Phase Portrait 

The quahtative behaviour of solutions is revealed by the phase portrait. See, for example, 
Figure[I], which is a phase portrait for certain values of the parameters. To find the horizontal 
and vertical isoclines, set, respectively, y = and i = in ([T]) to obtain the graphs 

y = H{x) := and y = V{x) :- 



1 + X 1 — r] + X 

Note that the EA corresponds to the vertical isocline and the QSSA corresponds to the 
horizontal isocline. (One may refer to y = H{x) as the quasi-steady-state manifold and 
y = V{x) as the rapid equilibrium manifold.) Observe that 

lim H{x) = 1 = lim V{x) and H{0) = = 1^(0). 

X— >oo x—^oo 

Since both H and V are strictly increasing and V{x) > H{x) for all x > 0, there is a narrow 
region between the isoclines: 

To := {{x,y) : x> 0, H{x) <y< V{x)} . 

Theorem 1. Let x(t) be a solution to ([T]). Then, there exists a t* > such that 

{x{t),yit))eToyt>t\ 
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Proof: First, we show that solutions can enter Tq but not leave it, that is, show that Fq is 
positively invariant. It follows from the differential equation that g • n < along both the 
vertical and horizontal isoclines, where n is the outward unit normal vector. Hence, Fq is 
positively invariant. 

Second, we establish that solutions outside Fq eventually enter Fq. Call y{x) the cor- 
responding one-dimensional solution. If y{x) is below the horizontal isocline H{x), then 
— < y'{x) < 0, and so y must intersect H for a lower value of x. Similarly, if y{x) is 
above the vertical isocline V{x), then — oo < y'{x) < —s~^, and so y must intersect V for a 
higher value of x. □ 

Behaviour of solutions near the origin, the only equilibrium point, is governed by the 
linearization matrix 



OX 



x=0 



-1 1-r] 

-1 -E-' 



(4) 



The eigenvalues of A are given by 



A± := , (5) 

which are real- valued and distinct. One can prove that 

A_ < -1 < A+ < 0, 

thus implying that the origin is asymptotically stable. Corresponding eigenvectors are 

Observe that v+ points into the positive quadrant while v_ does not. The slope of the 
eigenvector v_|_ at the origin is very important, and will be denoted by 

A+ + 1 

a := . 

I — f] 

Asymptotically, we have that 

a = l + er] + (e^) as e 0. 

Observe also that the slope of the slow manifold at the origin lies between the slope of the 
horizontal isocline and the slope of the vertical isocline. That is, 

1 < a < {l-ri)-\ 

The original, time- dependent differential equation ([T]) has linearization 

X = Ax, (6) 
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where the matrix A is as in (jlj) . Since the eigenvalues are real- valued and distinct, the initial 
value problem 

X = Ax, x(0) = xq 

has solution of the form 

x(t) = c_e^-*v_ + c+e^+*v+. 

We will assume, to avoid triviality, that xq 7^ 0. The coefficients c± can be determined in 
terms of left eigenvalues and eigenvectors. The left eigenvectors can be taken to be 

Using the orthogonality condition v_|_v=|i = we see that 

v±xo 
c± = ^ — • 

Proposition 2. Let y he a solution to ([3]) which lies inside Tq for x G (0,a), where a> 0. 
Then, 

lim y{x) = and lim y'{x) = a. 

x-^0+ x-s>0+ 

Proof: We should begin by emphasizing that solutions x(t) to ([1]) enter and forever remain 
in the interior of Tq. By hypothesis, H{x) < y{x) < V{x) for all x G (0,a). The Squeeze 
Theorem establishes the first limit since i^(0) = and ^(0) = 0. 

To establish the second result, observe that the function g, as in ([2]), is of class with 
g(0) = and the matrix A has strictly negative eigenvalues. It follows from Hartman's 
Theorem (see, for example, [21] p. 127), which is a stronger version of the Hartman-Grobman 
Theorem and applies even in cases of resonance, that the phase portrait of ([T]) behaves like 
the phase portrait of ([6]) diffeomorphically in a neighbourhood of the origin. Therefore, 
solutions to the nonlinear system have slope a as they approach the origin too. □ 

Remark 3. Since any solution x(t), except the trivial solution, eventually enters Tq and 
then approaches the origin asymptotically, we can now say definitively that the origin is 
globally asymptotically stable. 



3 The Isocline Structure 

The nature of the level curves, or isoclines, c = f{x,y{x)) reveals to us a surprising amount 
of insight into the behaviour of solutions to ([3]). Consider 

c=f{x,y{x)), (7) 

where / is as in the differential equation ([3]) and c G M. For a given x > and c G M, ([7j) is 
invertible and we can solve for y{x), yielding 

y(x) = F(x,c):=-^^, (8) 
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Figure 2: Graph of the function K{c) for arbitrary e > and 77 G (0, 1). 



for c —e 



, where 



K{c): 



l + £{l-r])c 
1 + ec 



Observe that y'{x) = f{x,y{x)) and y{x) = F{x,y'{x)) for solutions y of ([3]) for values of x 
for which the solution is defined. Throughout this paper, level curves of / will be denoted 
by w. If the slope associated with w is required, we specify this and write w{x) = F{x,c). 
For completeness, we will agree that F {x, —s^^) = 0. 

The function K, which is sketched in Figure [21 and the isoclines, which are sketched in 
Figure [3l have the following, easy-to-prove properties. 

Proposition 4. The isoclines and the function K satisfy the following. 

(a) The function K is strictly decreasing everywhere except at c = —e^^, where it has a 
vertical asymptote. 

(b) The function K satisfies 1 — rj < K{c) < 1 for < c < 00, which corresponds to the in- 
terior o/Fq. Furthermore, we have that K{0) = 1 corresponds to the horizontal isocline 
and linic^oo K{c) = 1 — rj corresponds to the vertical isocline. 

(c) The function K satisfies the important relation K{a) = a^^ . 

(d) Define the function u{c) := cK{c) for c > 0. The function u is strictly increasing, satis- 
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y 




Figure 3: Sketch of the isochne structure of ([3]) for arbitrary e > and t) G (0, 1). 



fits u{a) = 1, and 

.^ ^rN^/(0'l)' for ce{0,a) 
u{c) = cK{c) e <^ . 

1^(1, cx)), for ce{a,oo) 

See Figure\^ 

(e) Any of the isoclines w satisfy the differential equation 

w{w - 1) + xw' = 0. (9) 

Note that the isochnes are hyperbolas. There are two exceptional isoclines, namely 
w{x) = F (x, —e~^) = and w{x) = F (x, —£^"'^(1 — rj)^^) = 1. The vertical isocline, V, also 
is somewhat of an exceptional case. Approaching it from below, one encounters increasing c 
up to +00. After passing through V, the slopes increase from — oo. 

4 Existence and Uniqueness of the Slow Manifold 

We provide a brief review of fences and antifunnels, which form the backbone of the existence- 
uniqueness proof that follows. 

Phase spaces of differential equations often exhibit curious curves and regions known as 
fences, funnels, and antifunnels. The best source of information on funnels and antifunnels 
is [12], Chapters 1 and 4. 
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y 



y = u{c) 

y = {I- 7])c 
" c 

Figure 4: Graph of the function u{c) given in Proposition Hl^d) for arbitrary e > 
and 77 G (0,1). The function is strictly increasing for c>0 for every admissible e and 
rj. Furthermore, u'{0) = 1 and asymptotically u has slope 1 — 77 as c — i- 00. Incidently, 
u{c) — (1 — ri)c = e~^rj + O {c~^) as c — )■ 00. 



Definition 5. Let / = [a, h) be an interval (where a < b < 00) and consider the first-order 
differential equation y' = f{x, y) over /. Let a and /? be continuously-differentiable functions 
satisfying 

a'{x) < f{x,a{x)) and f{x,l3{x))<f3\x) (10) 

for all X G /. 

(a) The curves a and /3 satisfying (fTOj) are, respectively, a lower fence and an upper fence. 
If there is always a strict inequality in ( ITOj) . the fences are strong. Otherwise, the fences 
are weak. 

(b) If /3(x) < a{x) on /, then the set 

r := {ix,y) : X e I, (3{x) <y< a{x)} 
is called an antifunnel. The antifunnel is narrowing if 

lim \a{x) — (5{x) \ = 0. 

Theorem 6 (Antifunnel Theorem, [12] p. 31-33). Let T be an antifunnel with strong lower 
and upper fences a and (i, respectively, for the differential equation y' = f{x,y) over the 
interval I := [a,b) (where a < b < 00). Then, there exists a solution y{x) to the differential 
equation such that 

f3{x) < y{x) < a{x) for all x E I . 
If, in addition, T is narrowing and ^{x,y) > m F, then the solution y{x) is unique. 

We cannot use Fq as an antifunnel (in the sense of Definition [5]). The key to our proof 
is considering the isocline for slope a, the slope of the slow manifold at the origin. We will 
call this isocline a. Proposition Hl^c) and tell us that a is given by the simple expression 




The defining feature of a is cr = f{x,a{x)). Moreover, this function a has the remarkable 
property that a'(0) = a. That is, the isocline for slope a has slope a at the origin. Define 
the region 

Ti := {{x,y) : x > 0, H{x) <y< a{x)} , 
which is a subset of Fq because H{x) < a{x) < V{x) for all x > 0. 

Theorem 7. 

(a) There exists a unique slow manifold y = Ai{x) in Ti for the differential equation 

(b) The slow manifold y = A4.{x) is also the only solution that lies entirely inside Fq. 

Proof: 

(a) We will show that the Antifunnel Theorem can be applied to the interval [a, oo), where 
a > is arbitrary. First, we show that the curve y = a{x) is a. strong lower fence and 
the curve y = H{x) is a strong upper fence for the differential equation for x > 0. 
Now, the derivative of solutions along the concave-down curve y = a{x) is identically a. 
Thus, 

a'{x) < a'(0) = a = f{x, a{x)) V x > 0. 

Hence, by definition, y = a{x) is a strong lower fence for x > 0. To show that y = H{x) 
is a strong upper fence for x > 0, consider that 

f{x,H{x)) = < H'{x) Vx>0. 

Second, observe that the strong fences satisfy a{x) > H{x) for x > and 

lim \a{x) — H{x) \ = 0. 

X—>OD 

By definition, Fi is a narrowing antifunnel. 

Finally, a quick calculation shows that > in Fi. So, all the conditions for the 
Antifunnel Theorem (Theorem[6]) have been established. Therefore, there exists a unique 
solution y = Ai{x) to 1^ that lies entirely in Fi. 

(b) Obviously, any solution other than the slow manifold eventually leaves Fi. If the solution 
leaves Fi through the horizontal isocline it also leaves Fq, since both regions share the 
same lower boundary. If the solution leaves Fi through the a isocline, while in Fq the 
solution will have slopes in the range a < y' < oo and hence will eventually leave Fq 
since the upper boundary of Fq is bounded above by the line y = 1. 

□ 
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Remarks 8. 

(i) Theorem [7] shows that 

^ < M{x) < — V X > 0. 
1 + X a + X 

Thus, the necessity of the EA is diminished in the sense that a serves as a smaller 
upper bound on Ai. Furthermore, it follows from the isocline structure that A4{x) is 
strictly increasing, since solutions of the differential equation inside the antifunnel but 
not on the boundary have strictly positive slope. Note that this bound is especially 
tight when e is small, since 

a-^ = l-er] + (e^) as e ^ 0. 

(ii) Slow manifolds, like centre manifolds, are generally not unique and are defined locally. 
In our case, all solutions that have slope a at the origin are slow manifolds. However, 
we look at the global phase portrait and refer to the unique solution within Fi as the 
slow manifold. 



5 Concavity 

Let y be a solution to ([3]), which we assume is not the slow manifold because we will deal 
with that case later. Then, of course, y'{x) = f{x,y{x)) and so by the Chain Rule, 

y"{x) = p{x,y{x))h{x,y{x)), (12) 

where 

p{x, y) := e^^r] [—x + [1 — rj + x)yY'^ 

and 

Kx,y) ■.= y{y -I) + xf{x,y). (13) 

The function p{x^ y) is positive everywhere except along the vertical isocline, where it is 
undefined. The function h{x) := h{x,y{x)), the sign of which determines that of y"{x), has 
derivative 

h'{x) = 2y{x)y'{x) + xp{x)h{x), (14) 

where p{x) := p{x,y{x)). The concavity of all solutions in all regions of the non-negative 
quadrant can be deduced using this auxiliary function h. Table [1] summarizes what we will 
develop in this section. They are all suggested by the phase portrait in Figure [TJ 
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Region 


Concavity of Solutions 


0<y <M 


concave down 


M < y < a 


concave down, then inflection point, then concave up 


a <y <V 


concave up 


V <y <1 


concave down 


y>i 


concave up, then inflection point, then concave down 



Table 1: A summary of the concavity of solutions of ([2D in the non-negative quadrant. 



Remarks 9. 

(i) Let y he a solution to ([3]) and fix xq > 0. Define w{x) := F{x,y'{xo)) to be the iso- 
cline through {xo,y{xo)). By virtue of the isocline structure, y"{xo) > if and only 
if y'{xo) > w'{xo) and y"{xo) < if and only if y'^xo) < w'{xq). Indeed, from ^ and 

h{xo) = xo[y'{xo) - w'{xo)], (15) 

which confirms this fact. The similarity of the form of h{x) and the differential equation 
([H]) that the isoclines satisfy is not a coincidence. 

(ii) The function h cannot tell us anything about the concavity of solutions at x = 0, not 
even by taking a limit. 

Many of the following proofs will involve the following elementary lemma, so we single it 
out here. We omit the proof in the interest of space. 

Lemma 10. Let I be one of the intervals [a, b], (a, b), [a, b), and (a, b]. Suppose that (j) & C{I) 
is a function having at least one zero in I. 

(a) If I = {a,b] or I = [a,b], then the function (p has a right-most zero in I. Likewise, if 
I = [a, b) or I = [a, b], then the function (p has a left-most zero in I. 

(b) If (f) & C^{I) ^'^^ > for every zero of (p in I, then has exactly one zero in I. 

Proposition 11. Let y{x) be any solution to Q lying below the slow manifold, say with 
domain (0, a] and y{a) = 0. Then, y is concave down for all x G (0, a]. 

Proof: Let h be defined as in f|T3|) with respect to the solution y. There are two regions to 
consider, namely where y{x) > H{x) and where y{x) < II{x). It is clear from f|T2l) and f|T3l) 
that y"{x) < for y{x) < II{x), noting that < y{x) < 1 and y'{x) < 0. The solution y{x) 
crosses the horizontal isocline, say at x = X2 G (0,a). Here, h{x2) < using f[T^ . Suppose 
that the proposition is false and that there are one or more inflection points. Applying 
Lemma [TOl let xi G (0,X2) be the right-most zero of h. Now, from (fT^ . 

h'{xi) = 2y{xi)y'{xi) > 0. 



12 



Then, h is positive in a neighbourhood to the right of Xi. Since h{x2) < 0, by the Intermediate 
Value Theorem, h has a zero in (xi, X2) which contradicts the fact that Xi is the right-most 
zero. Therefore, there is no inflection point. □ 



Proposition 12. Let y be a solution to (E]) between a and V over {a,b). Then, y is concave 
up on (a, b). 

Proof: Fix xq G (a, b) and let c := y'{xo) and r := K{c). Let w{x) := F{x, c) be the isochne 
through {xo,y{xo)). With h defined as in ( !T3l) with respect to y, we have 

/i(xo) = xo[y'{xo) - ii;'(a;o)] = 

where we used the expression for h in (|T5l) . Since c > cr, applying Proposition H] we know 
rc > 1 which implies r > y/rcr^. Suppose, on the contrary, that y"{xo) < 0. Then, 

c(r + Xo)^ — r < Xq < V rc~^ — r < 0. 

This is a contradiction. □ 

Proposition 13. Let y be a solution to ([3]) lying between M. and a, with domain (0, a] and 
y{a) = a{a). Then, y has exactly one inflection point Xi G (0,a). Moreover, y is concave 
down on (0,xi) and concave up on (xi,a). 

Proof: We know y'{0) = a and y'{a) = a. Hence, by Rolle's Theorem, y has an inflection 
point Xi G (0, a). To prove uniqueness of the inflection point, let h be as in (fT3ll with respect 
to the solution y. Now, if x is a zero of h, then 

h'{x) = 2y{x)y'{x) > 0. 

By Lemma fTOl there is at most one zero of h. Moreover, since h'{xi) > 0, y is concave down 
on (0,xi) and concave up on (xi,a). 

□ 

Proposition 14. Let y be a solution to which lies above V and below 1, with domain 
[a, b) and lim^^f,- y{x) = V{b). Then, y is concave down for all x G [a, b). 

Proof: It is clear from the expression for h, (1131) . where h is defined with respect to the 
solution y, and the fact that < in that region. □ 

Proposition 15. Let y be a solution to ([3]) which lies above 1, with domain [0,a], where 
y{a) = 1. Then, there exists a unique inflection point xi G (0,a). Moreover, y is concave up 
over [0,xi) and concave down over (xi,a]. 



c(r + Xq) — r 
(r + xo)^ 
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Figure 5: The two thick curves are curves along which solutions of ([3]) have inflection points, 
for parameter values e = 0.6 and rj = 0.9. The thin curves are the horizontal, a, and vertical 
isoclines. 



Proof: Let h be defined as in ( fT3l) with respect to the solution y. Now, 

y'iO) = -e-\l-r^)-'=y'{a). 

By Rolle's Theorem, there exists xi G (0,a) such that y"{xi) = 0. The uniqueness of the 
infiection point follows from the fact that any zero x of h satisfies h'{x) < and an application 
of Lemma [TOl Moreover, since h'{xi) < 0, y is concave up on [0,xi) and concave down on 
{xi,a]. □ 

Remark 16. We now know that solutions can only have infiection points between Ai and a 
and above y = 1. There are, in fact, curves along which solutions have zero second derivative. 
To find them, one could, for example, use Maple to solve ^ {f{x,y{x))) = for y{x), the 
solutions unfortunately being rather long and messy. There are three solutions. One curve 
lies below the x-axis and is discarded. The other two curves are in the positive quadrant, 
one lying between Ai and a (and which is a lower fence actually), the other starting at (0, 1) 
and increasing with x. See Figure [51 



6 Behaviour of Solutions Near the Origin 

It was argued in pL9j and [22] , for example, that the slow manifold can be written as a Taylor 
series of the form J\4 (x) = CnX" at the origin. This is a traditional approach but we 

will show that this approach is not always valid. However, in the realm that is usually 
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considered for the Michaelis-Menten Mechanism, namely < e ^ 1, a very high number of 
terms of this Taylor series is correct. 

Intuitively, we know that A4 lies between the horizontal and vertical isoclines which both 
have limit zero as x — )■ 0+, and Ai shares the same direction as the slow eigenvector v+ at 
the origin. Hence, it must be that 

o"o = and o"i = a. (16a) 

By substituting the series into the differential equation, one can obtain all the coefficients 
recursively: 

Y2Zl [{n - k)an-k + {l-y]){n-k + l)an-k+i] cfk 

0"n = 



e ^ + (1 — vi){n + l)ai — n 



£ ^ + (1 — ?7)(n + l)ai — n 



(16b) 



Let y be any solution to ([3]) that lies inside Fq. Since no property of the slow manifold was 
used in constructing the above series which all other solutions do not possess, we can equally 
well write 

oo 

y{x) = Y,^^^'■ (17) 

Define 

.:^^^i±i±V^^S, (18) 

where we made use of the expression for A_ and A+, (jS]). 

Re-arranging (|T8l) . we see that t] can be written in terms of e and k as 

This can tell us when the parameter n takes on certain values. However, for a given e > 
and n G N\{1}, there may not be a corresponding rj G (0, 1) that gives k = n. It can be 
shown that 

^ = ^ + + C'(e) as £^0 

er] rj ^ ' 

and thus k — > oo as e — 0. That is, if e is very small, which is the case traditionally 
considered, k is very large. Many results that follow will involve k and so it is a good 
idea to keep this in mind. 

Observe that n > 1 and that we can choose values of the parameters e and r] to achieve 
any desired value of k we wish. The following is easy to prove. 

Proposition 17. Consider the constant k> 1 and the coefficient 

(a) There is resonance with the eigenvalues {A_,A+} if and only if n E N\{1}. (See, for 
example, [^J for a discussion of resonance.) 
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Figure 6: The graph of r]{e, k) for an arbitrary, fixed e > 0. The physically relevant values of 
K are k > 1 . Observe that > for < k < 1 and |^ < for «; > 1 with a global maximum 
at K = 1. Furthermore, there is an inflection point at k = 2 and — )■ as k — )■ oo. Observe 
that the maximum value satisfies > 1 for all e > and hence for any permissible value 

of e there are values of n (in a neighbourhood of /c = 1) which give inadmissible values of r}. 



(b) For the numbers k and a^, k, G (1,2) if and only if a2 > 0. Furthermore, k, > 2 if and 
only z/ (72 < 0. 

Proposition 18. Consider the constant k > 1 and the coefficients {cTij^Q- 

(a) If K ^N, then all o/{o"j}^o defined. 

(b) If K E N\ {1}, then {o-i}'^^^ are all defined but a^^ is not defined (and hence all subsequent 
cr„ are not defined). 

Proof: We know from f ll6ap that the coefficients ctq and di are always defined. Consider 
the expression (flBbl) . which gives the recursive descriptions of the coefficients. Solving 

+ (1 - ri){3 + l)ai - J = 0, J G {2, 3, . . .} 

gives 

From ( |T9ll . this is true if and only if k = j. Hence, if k G N\ {1}, then {ai}'l~Q are all defined 
but (Tk is not defined. If k ^ N\ {1}, then all the coefficients are defined. □ 

A classic method of finding an asymptotic expression for a solution to a one- dimensional 
differential equation y' = f{x, y) is the power series method. Here, one assumes a solution of 
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the form y{x) = Yl'^=o ^nX^i substitutes into the differential equation, and arrives at recursive 
relationships for the coefficients which are then solved. However, this is not always reliable. 

Theorem 19. Consider the system of ordinary differential equations 

x=Ax+h{x), x{{]) = xq, a;=(^)GM^ (20) 

where the matrix A is Hurwitz (asymptotically stable), the vector field b is analytic with 
||6(a;)|| = O (||a;||^) as \\x\\ 0, and \\xq\\ is sufficiently small. Let the eigenvalues be A_|_ 
and A_, where A_ < A+ < 0, and define the ratio k, := ^ > 1. Suppose that k (i.e. no 
resonance) and the eigenvector satisfies {v^)i ^ 0. If x(t) is a solution to (l20l) which ap- 
proaches the origin in the slow direction and (for simplicity) is strictly positive for sufficiently 
large t, then 

y{t) = ^anx{t)'' + Cx{t)'' + o{x{tf) as t ^ oo 

n=l 

for some constants {cr^jj^^^ (which are independent of initial condition) and C (which de- 
pends on the initial condition). 

Proof: In order to derive the necessary asymptotic expansion for y{t) in terms of x{t), we 
make use of the linearized problem 

z = Az, z(0) = zq. (21) 

To avoid the trivial solutions, which have nothing to offer us, we will assume that xq, zq 7^ 0. 
Let x(t) and z(t) be, respectively, the unique solutions to (!20|) and ( |2T|) . both of which tend 
to the origin as time tends to infinity. We will not consider the initial conditions xq and zq 
to be independent so that the solutions x(t) and z(t) can be related. Furthermore, we need 
both ||xo|| and ||zo|| to be small. 

The solution to the linear problem z(t) can be written in the explicit form 

z(it) = c+e^+*v+ + c_e^-*v_, 

where c+ > (since we assumed that solutions approach the origin from the right in the slow 
direction). 

We know that there is no resonance with the eigenvalues. Moreover, the eigenvalues are 
in the Poincare domain. Applying Poincare's Theorem (see, for example, |2] p. 190), there 
is a quadratic vector field q such that x = z + q(z). Hence, we can write (not uniquely if 
KeQ) 

x{t) ~ a^„e("^-+"^+)* = a„„e(™"+")^+* as t ^ 00 (22a) 

{m,n)(iS {m,,n)&S 

y{t)^ Y &™ne(™'-+"'+)* = Y as t ^ 00, (22b) 

{m,n)(^S (m,'n)(iS 
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where 

S := {{m,n) : m,n&'L, m,n>0, m + n > 1} . 

Let £ := [kJ . Then, the first £ + 1 most dominant terms in fl22|) are, in order of decreasing 
dominance. 

To see why this is the case, we make two observations. First, the fact that the listed 
exponentials are in decreasing order of dominance is obvious except maybe for the last two. 
Since k = ^ > £, we have A_ < £A+ < 0. Second, there cannot be any other exponentials 
of the form e("^^+ in between those listed. 

For our purposes, we need only the first £ + 1 terms of ([22]) and hence we write 

e 

x{t) = ame"'^+^ + a^+ie^-* + o (e^-*) as t ^ oo (23a) 

m=l 

= J2 bme""^+' + be+ie^~' + o (e^-*) as t ^ oo, (23b) 



m=l 



where ai ^ 0. The coefficients can be related using the differential equation. To write y{t) 
in terms of x{t), we will successively eliminate the exponentials. Manipulating f l23ap and 

y{t) - a^x{t) = J2 + bfl.e'^' + o (e^"*) , (24) 

m=2 

where ai := — and bm := — (Jiam- To go further, observe that we can write powers of 
x(t) as 

e 

x(t)" = <e"^+* + J2 Cmne"'^+' + o {e^-') as t ^ oo, 

m=n+l 

where n G {2, . . . ,£}. Solving for the most dominant exponential. 



^a;(t)"- J2 (^) e™^+* + o (e^-*) as t ^ oo. (25) 



m=n+l 



With f l25|) . we can successively eliminate the exponentials of flMj) — each time introducing 
other exponential terms but none of order already eliminated — until we are left with an 
expression of the form 

e 

y(t) - ^cT„x(t)"^ = 6f+/^e^-* + o(e^-*) as t ^ oo. (26) 

m=l 

Since x{t)'^ = ale^-^ + o (e^-*) and hence 

e^-* = j^xity + (e^-*) as t -> oo, 
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we can write ( 126|) as 



y{t) - amxit)"" = Cx{tY + o (e^-*) as t ^ oo, 

m=l 

where 

^ bi+i — cTiae+i 

The desired conclusion follows. □ 

Remark 20. The coefficients are calculated using the power series method. That 

is, one assumes that the solution to the one-dimensional version of (|2Ui) is y{x) = Yl'^=o ^nX'^ 
(where ctq = by necessity). The purpose of the theorem is to tell us how many of the 
resulting terms apply to all solutions. 

Corollary 21. There exists a solution to (12U]) such that 

oo 

~ (T„x(t)'^ as t — !• oo 

n=l 

for some constants {cr„}^]^. 

Proof: Choose the initial condition zq so that it is parallel to the slow eigenvector v_|_. Then, 

z(t) = c+e^+*v+ 

and hence we can write 

oo oo 

x{t) ~ ^ 0^6""^+* and y{t) ~ ^ &me™^+* as t ^ oo. 

m=l m=l 

Any positive integer power of x{t) will be a series of the same form as x{t) and y{t). Suc- 
cessively eliminating exponents, just like in the proof of the theorem, gives us our desired 
conclusion. □ 

Now, we apply this general result to the Michaelis-Menten Mechanism. 

Lemma 22. Let y be a solution to ([3]) lying inside Tq and let the ratio of the eigenvalues be 
K > 1. Suppose that k (i.e. there is no resonance). Then, 

= ^cr„x" + Cx'^ + o(a;'') as x -)■ 0+, 



n=l 



where {c"n}i=i 0''i"e as in (fT71) and C is some constant that distinguishes the solution y{ 
from other such solutions. 
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x] 



Proof: This proof is a simple application of Theorem [191 The original, time-dependent 
differential equation ([T]) can be written 

i = Ax + b(x), x(0)=xo (27) 

where the analytic vector field b is given by 

b(x) := xy ( ) . 

Observe that 

||b(x)|| = C (||xf) as ||x||^0. 
This is fairly obvious but can be shown directly: 

< (x - yy = + - 2xy = W^f - ^^\\hix.)\\ 

and hence 

||b(x)|| < vipEiixf . 

The linearized problem is 

z = Az, z(0) = zo, (28) 

where the matrix A = (^_\ ^^-i ) was given in (jlj). Again, we will assume that xo,zo 7^ 
and, in particular, lie in the positive quadrant. Let x(t) and z{t) be, respectively, the unique 
solutions to (127]) and fl28l) . The initial conditions xq and zq are not independent so that the 
solutions x(t) and z{t) can be related. Furthermore, we need both ||xo|| and ||zo|| to be small. 
The solution to the linear problem z(t), as we have seen earlier, can be written explicitly as 

z(t) = c+e^+*v+ + c_e^-V_, c± := 

with c+ > and the sign of c_ depending on which side of v+ the initial point zq lies. 
Finally applying Theorem [191 after dropping the time dependence we can say that 



y{x} 



(T„a;" + Cx"" + (x") as x 0^ 



n=l 



for some constants {S'n}l=i and C, where i := [k\ . By uniqueness, we have a„ = cr„ for all 
ne{l,. ..,£}. □ 

Remark 23. For k G (1,2), we can manipulate the given constants to get 

^ ^ C-(A- - A+) 
cl{l - riY 
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7 All or Most Solutions Must Enter the Antifunnel 



We now investigate conditions under which solutions enter Ti. 

Theorem 24. Let x{t) be a solution to ([T]) and suppose there is no resonance, i.e. ^ N. 

(a) If K> 2, then there exists a t* > such that 

(b) If K, < 2, then there exist solutions x{t) which do not enter Ti. Moreover, solutions 
that do not enter Ti must enter Fq through the vertical isocline V to the left of the line 
y = ax. 

Proof: We begin by noting that if a solution x(t) enters Fi, it forever remains in Fi. 
This is because g • n < along a and H, where n is the unit normal vector. Let y be the 
corresponding one- dimensional solution. 

(a) Applying Lemma [22|, we can write 

y{x) = ax + a2X^ + o(x^) as x — )• 0^. 
Furthermore, since k > 2 we have (72 < 0. Thus, 

lim y"{x) = 2a2 < 0. 

x-s>0+ 

Since solutions are concave down only when they lie below the isocline a, x(t) eventually 
enters the Fi antifunnel. 

(b) From Lemma we have 

y{x) = ax + Cx" + a (x*^) as z -> 0+ 

for some constant C. It follows that there are some solutions to that are concave up 
at the origin — the ones for which C > — and curve away from a and exit Fq through 
the vertical isocline. Moreover, since ^2 > 0, by virtue of Corollary |2T] it follows that 
there is a solution with a Taylor series at the origin that does not enter Fi from above. 
See Figure [3 

We know already that x(t) eventually enters Fq. If x(t) enters Fq through the horizontal 
isocline, it also enters Fi. Denote by {x* , y*) the point of intersection of the line y = ax 
and the vertical isocline y = V{x) (that is, y* = ax* = V{x*)). Assume that x(t) enters 
Fq through the vertical isocline to the right of {x*,y*). We claim that x(t) also enters 
Fi. Suppose that y intersects the vertical isocline at x = xi. Observe that xi > x*, by 
assumption, and y{xi) < axi. By the Mean Value Theorem, there is a xq G (0, xi) such 
that 

0<y [xo) = — = < = a. 

Xi — (J Xi Xi 

By Virtue of the isocline structure, H{xq) < y{xo) < a{xo) and therefore x()f:) eventually 
enters Fi. 

□ 
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" y = ax 




• X 



Figure 7: Sketch demonstrating Theorem [2^b). 

8 Properties of the Slow Manifold 

Finally, we present some important properties of the slow manifold. 
Proposition 25. The slow manifold y = Ai{x) is concave down for all x > 0. 

Proof: Construct a sequence of functions {yn}'^=N follows. Fix xq > and let be the 
solution to ([3]) such that 

yn{xo) = M{xo) - 

The number is taken large enough so that yNi^o) > 0. Thus, by Proposition [TTl 2/^(xo) < 
for all n. Now, 

y'nixo) = pixo, yn{xo))h{xo, yn{xo)) (29) 

and 

M"ixo) = pixo,Mixo))hixo,Mixo)). (30) 

By construction, 

lim yn{xo) = M{xo). 
Letting n — )■ oo in fl29|) and applying (|30|) . continuity tells us 

lim y'^{xo) = M"{xo). 

n— ^oo 

Since y'nixo) < for all n, Ai"{xQ) < 0. Since Xq was arbitrary, Ai"{x) < for all a: > 0. 



22 



We now establish a strict inequality. Suppose that Ai"{x*) = for x* > 0. If /i is as 
in (IT^ with respect to the solution Ai, we have h'{x*) > 0. This contradicts the fact that 
h{x) < for all x > 0. □ 



Proposition 26. The slow manifold y = Ai{x) satisfies, for all x > 0, 

< H{x) < M{x) < a{x) < 1. 

Furthermore, 

lim M{x) = and hm M{x) = 1. 

Proof: The first result follows from Theorem [7] and the definition of a, (fTTj) . To establish 
the other two results, note that the functions H, a, and A4 are all continuous for x > 0. 
Since 

lim H{x) = lim a{x) = and lim H{x) = lim a{x) = 1, 
the results follow from the Squeeze Theorem. □ 

Proposition 27. The slope of the slow manifold y = M.{x) satisfies, for x > 0, 

< M'{x) < a. 

Furthermore, 

lim M'{x) = a and lim M'{x) = 0. 

x->-0+ x-s>oo 

Proof: The first result follows from the fact that the slow manifold lies within Fi, which 
consists of nested isoclines of slopes varying from to cr. 

To prove the second result, observe that the direction of the slow manifold at the origin 
must correspond to the slow eigenvector at the origin, which has slope a. 

To prove the third result, we note that, since Ai is strictly increasing and concave down, 
there is a c G [0, cr) such that 

lim Ai'{x) = c. 

x—^oo 

Suppose, on the contrary, that c > 0. By the Fundamental Theorem of Calculus, 

/•X PX 

M.{x) = / Ai'{u)du> / cdu = cx. 
Jo Jo 

However, for sufficiently large x, cx > a{x), a contradiction. □ 
The slow manifold has also been approximated, for large x, in the asymptotic series [22] 

Ai{x) ~ po + Pix~^ + P2X~^ + ■ ■ ■ as X — )■ oo. (31) 
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The coefficients can be obtained by substituting the series into the differential equation and 
are given recursively by 

Po = 1, Pi = -1, P2 = 1, 

n-2 

Pn = -pn-1 + ^ X] ipi[pn-i-i + (1 " v) Pn-i-2] for u > 2. (32) 

i=l 

Observe that all the coefficients are polynomials in e and r]. As we will establish now, the 
series (l3T|l is fully correct. 

Proposition 28. For large x, the slow manifold satisfies 

oo 
n=0 

where {p^j^Q are as in (!32|) . For small x, in the case of no resonance (i.e. k ^N), the slow 
manifold satisfies 

M{x) = '^aiX^ + Cx" + o{x'^) as a; 0+ 

i=l 

for some constant C . 

Proof: The second conclusion follows from Lemma [221 To prove the ffist conclusion, observe 
that for any c > there exists a x* > such that 

H{x) < M{x) < F{x,c) 

for all X > X*. This is because Ai is concave down and lim^_^oo = 0. Hence, 

1 - x~^ + O [x-"^) < M{x) < 1 - K{c)x-^ + O [x''^) as x ^ oo. 

Since c > was arbitrary and K{c) — )■ 1 as c — )• 0, it follows that 

x^oo — X 

Hence, 

Ai{x) = 1 — x~^ + (x^^) as X — >■ cxo. 

Unfortunately, this is as much information that we can extract using the isoclines. To obtain 
the remaining terms of the asymptotic series, we will use the Centre Manifold Theorem (see, 
for example, [4j). 

Under the change of variables 

X:=x-\ Y ■.= y- {1-x-^) , (33) 
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we arrive at the system 



X = -XV {X-\l-X + Y) 

Y = ~X^g, {X-\ l-X + Y)+g2 {X-\ 1-X + Y), (34) 

where gi and g2 are as in ([2]). The system is not polynomial but there is no harm, 
because the resulting one- dimensional differential equation will be the same, in considering 
the system 

X = -X'gi {X-\l-X + Y) 

Y = -X'g, {X-\ 1-X + Y)+Xg2 {X-\ l-X + Y), (35) 

which is polynomial. Expanding, we get the expressions 

X^gi {X-\ 1-X + Y) = -X^ [riX-Y+{l-ri)X{X- Y)] 
Xg2 {X-\ 1-X + Y) =£-^ {X^ -XY -Y) . 

The system (!35ll . as we see, is a bit more messy than the original system ([3]). The eigenvalues 
of the matrix for the linear part of the new system ( |35l) . which is diagonal by construction, 
are and — e:~^. We know from centre manifold theory that there is a centre manifold which, 
we claim, must be the slow manifold. 

Observe that the "K-axis is invariant. Moreover, the fixed point {X, Y) = (0, 0) is a saddle 
node (or a degenerate saddle). The physically relevant portion of the phase portrait, namely 
X > and y > — 1, consists of two hyperbolic sectors, one with the positive Y-axis and 
the centre manifold as boundaries and the other with the negative y-axis and the centre 
manifold as boundaries. See Figure El This can be shown using techniques in §9.21 of [Ij 
(in particular Theorem 65). This also is a consequence of the phase portrait of the original 
system ([T]). It follows that the centre manifold is indeed the slow manifold. 

By the Centre Manifold Theorem, the slow manifold (in the new coordinates) can be 
written 

oo 

M{X)^J2 Pi^' as X ^ 0+, 

i=2 

for some coefficients {pj}^2- '^^^ ^^^^ terms are given by p2 = 1 and 'p^ = ei] — 1. 
Reverting back to the original coordinates, 

oo 

Ai{x) 1 — x~'^ + ^^piX"'' as a; — 7- oo. 

1=2 

Observing that the coefficients in (132|) are generated uniquely from the differential equation, 
the conclusion follows. □ 
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Figure 8: A phase portrait for ( l35l) for e = 5.0 and 77 = 0.8. 



Remarks 29. 

(i) We use the ad hoc transformation fl33|) because it is inspired by the series for Ai 
which we wish to obtain and it also resuhs in a system which is in the canonical form 
of the Centre Manifold Theorem. Others, for example, P, |T0], have used Poincare 
compactification to study the behaviour of Ai at infinity and found that the fixed 
point is a degenerate saddle. 

(ii) The Centre Manifold Theorem can be applied at the origin as well. See, for example, 
[1] pages 8-10. However, this result gives a smooth solution for small e only. This is 
because in order to apply the Centre Manifold Theorem, the differential equation 6 = 
is appended to the system ([2]) which gives the zero eigenvalue. A centre manifold exists 
in a neighbourhood of {x, y, e) = (0, 0, 0). 

Proposition 30. In the case of no resonance, i.e. k ^ N, the second derivative of the slow 
manifold satisfies 



Proof: The proof involves an easy application of Lemma [22] and the fact that the slow 




manifold is concave down at the origin. 



□ 



26 



Remark 31. The quasi-steady-state approximation has traditionally been used to approxi- 
mate the long-term behaviour (in time) of solutions to ([T]). Is this justified? Recall that 



H(x) = and a(x) = — ; . 

It follows that the QSSA is good when a ~ 1. Recall also that a = 1 + O (e) as e — t- 0. 
Hence, the QSSA is a good approximation when e is small. However, the function a{x) has 
slope a'(0) = o" at the origin and so is a good approximation for solutions near the origin for 
any e. Furthermore, since k = [sr])'^ + as £ — )■ 0, a large number of Taylor coefficients 
are correct in the asymptotic expansion at the origin for the slow manifold if e is small. 



9 Open Questions 

In the analysis of the behaviour at the origin, we assume non-resonance of the eigenvalues 
of the linearization at the origin. The resonance cases still need to be investigated. 

S. Fraser and M. Roussel[Hl UHl |22], [231 El] have introduced and investigated an iter- 
ation scheme to approximate the slow manifold. Specifically, the iterates are defined by 
i/n+i '■= F{x, y^), where F is as in ([8]). A definitive proof of convergence of the scheme which 
is valid for all values of the parameters e and i] has not yet been given. The scheme can 
diverge for certain values of the parameters and certain choices of initial iterate and converge 
for others. Convergence has been examined, for example, in [13]. In this particular paper, 
the Fraser iterates and perturbation series for JH in e were compared. Specifically, if {yn}'^=o 
are the Fraser iterates with initial iterate Uq := H and 

oo 

M{x) = Y,Mm{^)e"' 

m=0 

is the perturbation series for Ai, then for each n 

n 

yn{x) = J2-Mm{x)e'^ + 0(6''-^') as £->0. 

m=0 
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